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This  paper  is  concerned  vlth  an  empirical  evaluation  of  the  relative 
performance  of  several  selected  edge  detectors  for  the  detection  of  edge  struc- 
ture in  noisy  digitized  images.  In  particular,  ve  consider  both  the  Hueckel 
operator  and  a recently  introduced  class  of  edge  detectors  Implemented  as  two- 
dimensional  infinite  Inpulse  response  (HR)  or  recursive  digital  filters.  The 
latter  were  originally  developed  on  the  basis  of  least  mean-square  Wiener 
spatial  filtering;  concepts  for  an  assumed  stochastic  model  of  edge  structure 
in  typical  imagery  data.  These  two  schemes  are  compared  with  respect  to  per- 
formance, computational  complexity  and  suitability  for  specific  applications. 
Problems  in  making  these  comparisons  are  discussed.  The  Implementations  of 
both  operators  used  for  this  comparison  are  described  in  detail. 
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I.  Introduction: 

A recurring  problem  in  pattern  recognition  applications  is  the  detection 
of  edges  in  digitized  images.  Here  ve  utilize  the  intuitive  notion  of  an  edge 
element  as  a picture  element  (pixel)  lying  on  the. boundary  between  two  objects 
or  regions  of  widely  different  gray  levels . The  requirements  of  a good  edge 
detector  in  this  context  are  that  it  respond  only  to  true  edge  structure  and 
be  relatively  insensitive  to  noise  and/or  spurious  detail.  Additionally,  ve 
would  expect  that  the  edge  detector  possess  a computationally  efficient  reali- 
zation. Unfortunately,  very  little  information  is  available  concerning  the 
relative  performance  and/ or  computational  complexities  of  various  edge  detectors 
operating  on  vide  classes  of  imagery  data.  A possible  exception  is  the  work  by 
Bullock  [l]  although  even  here  a rather  restricted  class  of  images  were  consider- 
ed and  no  consideration  was  given  to  questions  of  computational  complexity. 

This  paper  is  concerned  then  with  an  empirical  evaluation  of  the  relative 
performance  of  several  selected  edge  detectors  for  the  detection  of  edge 
structure  in  noisy  or  textured  images.  In  particular,  ve  consider  both  the  Hueckel 
operator  [2], [3]  and  a recently  Introduced  class  of  edge  detectors [b]  implemented 
as  two-dimensional  (2-D)  infinite  impulse  response  (IIR)  or  recursive  digital 
filters.  The  latter  were  originally  developed  on  the  basis  of  least  mean-square 
Wiener  spatial  filtering  concepts  for  an  assumed  stochastic  model  of  edge  struc- 
ture in  typical  Imagery  data.  These  two  schemes  are  compared  with  respect  to 
performance,  computational  complexity  and  suitability  for  specific  applications. 

For  comparison  purposes  we  also  consider  the  U-point  Laplacian  operator.  The 
specific  implementations  of  both  the  Hueckel  and  Wiener  operators  employed  in 
this  study  are  described  in  detail. 
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Relative  evaluation  of  edge  detector  performance  is  of  necessity  some- 
what heuristic.  The  particular  application  generally  dictates  the  seriousness 
of  missed,  spurious,  or  excessively  thick  edges.  Indeed,  for  most  real-world 
images,  the  differentiation  between  true  edge  structure  and  spurious  detail  is 
somewhat  arbitrary.  Nevertheless,  it  is  expedient  to  rely  on  visual  evaluation 
criterion, however  subjective  this  may  be,  to  assess  the  accuracy  with  which  the 
edge  detector  output  represents  the  essential  "edge  information”.  This  is  the 
approach  taken  here. 

In  assessing  the  computational  burden  associated  with  various  edge  detectors, 
it  is  again  difficult  to  define  meaningful  quantitative  measures.  The  cost  of 
an  operation  depends  upon  the  precision  required  as  well  as  the  type  of  opera- 
tion (additlon/subtraction  vs.  multiplication/division).  The  relative  costs 
of  the  operations  varies  from  machine  to  machine  for  programs  written  at  the 
assembly  language  level  and  is  even  less  determinate  when  a higher  level  language 
is  used.  Also  the  inherent  structual  complexity  of  the  algorithm  must  be  con- 
sidered as  a simple  algorithm  lends  itself  more  readily  to  implementation  in  a 
special  hardware  processor.  To  the  extent  possible,  these  factors  have  been 
taken  into  consideration  in  the  assessment  of  relative  computational  complexity. 

Finally,  all  experimental  work  was  conducted  on  the  RPI  image  processing 
facility.  Basically,  this  is  a versatile  hardware/software  system  for  the  digit- 
isation, processing  and  display  of  image  data.  The  system  is  configured  around 
a Varian  620/1  minicomputer  with  32K  core  storage,  and  includes  a 5 megabyte 


disk  for  storage  of  image  files.  Additional  details  of  this  system  are  describ- 
ed in  [5].  The  software  support  is  coded  in  assembly  language  to  minimize  pro- 
cessing times  so  that  the  facility  can  be  used  interactively.  In  the  work  des- 
cribed here,  the  longest  processing  time  for  a single  operation  on  a 256  x 256 


raster  image  is  on  the  order,  of  a minute.  The  system  includes  locally  written 
routines  that  provide  efficient  I/O  for  fixed  length  random  access  files  (e.g.. 
Images ) . 

II.  Preliminaries ; 

In  choosing  a criterion  for  comparing  the  relative  performance  of  these  oper- 
ators, consideration  must  be  given  to  the  intended  use  of  the  output.  An  edge 
operator  which  responds  to  the  boundary  between  two  objects  by  indicating  many 
small  line  segments  which  are  close  together  and  parallel  but  not  necessarily 
connected  would  not  present  a serious  problem  to  a higher  level  processor  which 
was  capable  of  Interpreting  each  line  segment  in  context  with  surrounding  segments. 
To  a contour  tracing  algorithm,  however,  the  excess  of  marked  edge  points  plus 
the  large  number  of  alternate  paths  represent  a computational  burden.  We  are 
Implicitly  assuming  here  that  the  contour  tracing  algorithm  is  implemented  as 
some  form  of  heuristic  tree  searching  algorithm  such  as  described  in  16].  In 
addition,  if  a constraint  of  finite  storage  is  imposed,  the  multiplicity  of 

f 

potential  paths  so  created  makes  it  more  likely  that  the  algorithm  will  discard 
all  correct  paths  and  fail  to  extract  a complete  contour.  The  Hueckel  operator 
is  an  example  of  such  a processor;  its  results  are  based  on  data  in  a local  area 
and  need  further  Interpretation  in  order  to  provide  conclusions  about  the  content 
of  the  image.  The  Wiener  operator,  on  the  other  hand,  makes  decisions  about  the 
presence  of  an  edge  on  a point-by-point  basis  but  utilizes  observations 
over  the  entire  image  field.  It  has  been  designed  so  that  it  tends  to  divide 
the  image  into  contiguous  disjoint  regions  by  marking  as  edge  elements  continuous 
strings  of  points.  This  feature  is  especially  valuable  to  a contour  tracing 
algorithm  of  the  tree  searching  variety. 
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Obtaining  meaningful  measurements  of  the  computational  burden  associated 
with  an  algorithm  requires  careful  consideration  of  the  types  of  operations 
performed  and  the  precision  required  for  each  and  is  further  complicated  if 
conditional  branches  are  present.  The  programs  described  here  use  two  levels 
of  precision:  16  bit  and  31  bit  fixed  point  arithmetic.  The  following  table, 

although  somewhat  arbitrary,  has  been  used  to  provide  a measure  of  relative 
computational  costs: 


Operation 

Cost 

16 

bit 

addition  or  subtraction 

1 operation 

16 

bit 

multiply 

3 operations 

31 

bit 

add-subtract 

2 operations 

31 

bit 

multiply 

12  operations 

31 

bit 

square  root 

20  operations 

TABLE  1 

Assignment  of  Relative  Computational  Costs 


In  making  these  comparisons  it  should  be  kept  in  mind  that  a recursive 
algorithm  is  generally  easier  to  program  optimally  and  furthermore  is  more 
readily  implemented  as  a special  hardware  processor.  The  Hueckel  operator  is 
not  of  this  category.  When  described  in  a high  level  language  it  is  deceptively 
simple;  when  expressed  in  terms  of  machine  level  code  (i.e.,  assembly  language) 
it  becomes  very  lengthy  and  quite  complicated.  The  Wiener  operator,  on  the 
other  hand,  consists  of  a 2-D  recursion  relation  which  is  readily  written  in 
assembly  language  and  for  which  a special  processor  could  easily  be  built. 

The  method  of  uaeage  of  the  operator  must  also  be  taken  into  consideration 
as  situations  may  occur  where  information  is  needed  about  only  a part  of  the 
image.  The  Hueckel  operator  is  readily  applied  to  subsections  of  the  image 
since  each  application  uses  as  input  only  data  in  a finite  local  area.  Since 
the  Wiener  Operator  is  recursive  it  must  be  applied  to  the  entire  image  in 
order  to  obtain  information  about  any  part  of  it. 
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Ill . Wiener  Operator  Description: 


The  Wiener  edge  detector  is  derived  by  formulating  the  problem  of  edge 
detection  in  the  presence  of  noise  as  one  of  least  mean-square  error  spatial 
filtering  where  the  desired  operation  is  the  Laplacian.  This  requires  specif- 
ication of  a stochastic  model  for  images  possessing  inherent  edge  structure.  A 
useful  model  for  this  purpose  has  been  described  previously  in  [4]. . According 
to  this  model,  the  plane  is  partitioned  into  disjoint  rectangular  regions,  the 
sides  of  which  are  determined  by  a pair  of  mutually  independent  Poisson  point 
processes  evolving  along  each  of  the  orthogonal  coordinate  axes.  Gray  levels 
are  assigned  to  these  elementary  rectangles  according  to  a 2-D  discrete  Gauss- 
Markov  random  field.  Two  descriptive  parameters  are  available:  X represents  the 
edge  density  in  events  per  unit  dlstance^while  p represents  the  correlation 
coefficient  of  gray  levels  in  contiguous  regions. 

The  optimum  Wiener  filter  is  completely  specified  in  terms  of  these  two 
parameters  in  addition  to  the  quantity  £ which  represents  the  SNR  of  the  edge 
structure  process  vis-A-vis  the  white  background  noise  and/or  Bpurious  detail. 

As  shown  in  [4]  the  resulting  Wiener  filter  is  isotropic  with  system  transfer 


function 


H0(r) 




[r’.gll-piyi  AMl-pW 
6(1-P)X£ 


; r > 0 


where  r*  • and  u>^,  i»l,2,  represents  a spatial  frequency  component. 

The  Wiener  filter  is  implemented  in  the  discrete  domain  as  a 2-D  IIR  digital 
filter  (cf.  (7]«  [8]).  Briefly,  the  system  transfer  function  of  this  filter  is 


described  by 


H<*1**2>  “ *1^*2  *1**2^ 


t Actually  in  what  follows  we  assume  X is  measured  in  events/pixel. 
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where 


H1(z1»z2)  * H0(z1»Z2)+Z1H0(z11’Z2)+Z2H0(z1’Z21)+Z1  'W  Zi‘ ’*2*  > (3) 


and  is  a simple  first-order  section  of  the  form 


H0(Z1»Z2)“A 


1 - l^fbjj+lHz"  +z"  X+bjjz"  z" 

1 + aio<zI1+^1)+aiiVZ21 


(U) 


which  can  be  implemented  by  the  recursion  relation 

A computer  program  has  been  written  for  determining  the  three  coefficients 
a , aM  and  bM  and  the  gain  A according  to  an  iterative  gradient  procedure 
to  result  in  a frequency  response  for  H (z^,Zg)  which  provides  a least  mean- 
square  approximation  to  the  desired  response  HQ(r).  The  details  of  this  pro- 
gram are  described  in  [9].  In  Table  2 we  summarize  the  results  of  this  proced- 
ure for  selected  values  of  p,  X,  and  £. 

Finally,  one  further  stage  of  processing  1b  performed  before  declaring  an 
edge  present.  The  output  of  the  Wiener  filter  is  examined  in  3 x 3 arrays 
centered  upon  the  point  under  examination  at  position  (i,j).  Each  of  the  four 
syanetric&lly  opposite  pairs  of  points  are  examined  and  the  pair  corresponding 
to  the  largest  magnitude  of  the  quantities 

Ai  " ~ Xi+i,j+i 

“2  - ^ (6) 

■ f’W’W.j’  • 1Mk 

is  selected  end  a tentative  direction  code  is  assigned  to  the  point  on  the  basis 

of  the  pair  of  points  chosen  and  the  sign  of  their  difference.  Finally  the  above 

• ' ■'  v ■ :* " 

direction  code  aad  magnitude  are  accepted  as  marking  an  edge  if  for  specified 
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bu 

-0.825 6 

-0.8921* 

-0.8252 

-0.5736 

-0.8261* 

-0.8256 

■ 

X-0.0125 

*10 

-0.2ll»9 

-0.2565 

-0.3217 

-0.1*651* 

-0 . 3602 

-0.381*0 

-0.1*815 

*11 

0.0098 

-0.2189 

-0.1579 

0.01*1*9 

-0.151*0 

-0.11*61 

0.6102 

A 

0.115 

0.053 

0.038 

0.025 

0.022 

0.015  ’ 

0 . 008£ 

bu 

-0.8258 

-0.9396 

-0.81*60 

-0.8228 

-0.8661* 

-0.8120 

-0.826c 

-0.211*9 

-0.1989 

-0.281*9 

-0.3106 

-0.3075 

-0.3686 

-0 . 3782 

X-0.025C 

xo 

a 

11 

0.0098 

-0.21*90 

-0.1687 

-0.1761 

-0.201*0 

-0.1350 

-0.1531* 

A 

0.115 

0.07 

0.053 

0.01* 

0.031* 

0.023 

0.015 

*11 

-0.8256 

-0.921*1* 

-0.7980 

-0.81*80 

-0.8896 

-0.8228 

-0.7636 

-0.211*9 

-0.191*1* 

-0.2925 

-0.2821 

-0.2665 

-0.3326 

-0.3929 

X-0.050 

iU 

a 

11 

0.0098 

-0.1929 

-0.0881* 

-O.1696 

-0.2213 

-0.1575 

-0.0895 

A 

0.115 

0.081* 

O.O67 

0.051* 

0.01*8 

■ 

. 

0.031* 

0.023 

n 


Notes:  1.  For  filter  coefficients  are  independent  of  1 and  p. 

2.  Additional  results  are  available  from  the  authors. 


0 


TABLE  2 


parameters  F and  T the  following  inequality  is  satisfied 


my  A.  and  the  summation  is  over  the  3x3  array  centered  on  loca- 
l<i<l»  1 

The  computations  required  at  each  pixel  are  summarized  in  Table  3. 


16  bit  nmltiplication 


16  bit  addition 


TABLE  3 


Computations/Pixel  for  Wiener  Operator 


According  to  Table  1,  the  Wiener  Operator  requires  an  equivalent  of  105  operations 


per  pixel.  The  program  listing  in  Appendix  II  is  included  as  a description  of  this 


IV.  Hueckel  Operator  Description 


The  Hueckel  operator  has  been  described  in  [3]  and  [U],  Briefly,  the  oper- 
ator tries  to  fit  an  ideal  edge  element  to  the  image  data  in  the  domain  (a  circu- 
lar disk)  such  that  a minimum  mean-Bquared  error  criterion  is  satisfied.  An 


initial  step  in  this  procedure  is  the  least  mew-square  projection  onto  a set 
of  eight  basiB  functions  H. (x,y),  i*l,2,...,8  or  their  discretized  equivalents 


A computationally  efficient  algorithm  has  been  described  for  fitting  an  ideal 


edge  element  to  the  data  on  the  basis  of  this  projection 


The  Hueckel  operator  as  used  here  is  a hybrid  between  an  implementation 
employed  by  Bullock  [l]  and  that  described  in  [3].  The  resulting  algorithm  is 
described  below.  It  was  implemented  in  assembly  language  on  a Varian  620/1 
using  16  bit  arithmetic  to  compute  the  projections  A.,  i»l,2,...,8  onto  the 


corresponding  basis  functions.  Due  mainly  to  dynamic  range  problems  it  was 


avmni'ci 


found  necessary  to  employ  31  bit  fixed  point  arithmetric  in  fitting  an  ideal 


edge  element  to  the  projected  data.  The  calculation  of  the  projections  i*l, 
2,..., 8 takes  advantage  of  all  symmetries  present  in  the  basis  functions.  This 
part  of  the  operator  executes  in  approximately  2 5 seconds  for  a 256  v 256  image. 

In  computing  discretized  versions  of  the  basis  functions  the  circular  disk 
is  assumed  partitioned  into  69  square  regions  as  indicated  in  Fig.  1.  The  dis- 
cretized basis  functions  ( J ) , i*l,2,...,8,  J*l,2,...,69  were  determined  by 

numerical  integration  of  the  corresponding  H^Xjy)  over  each  of  the  69  square 
regions.  That  is,  for  an  interior  square  we  define 

H ( J ) = f [ H (x,y) dxdy  ; i-1,2 8 

1 J J 1 J=l,2, . . . ,69  (8) 

SJ 

where  Sj  is  the  square  region  of  integration  appropriate  for  the  J'th  square. 
Basis  functions  corresponding  to  squares  along  the  periphery  of  the  disk  have 
been  modified  to  include  edge  effects.  In  particular,  sections  of  the  disk 
not  in  any  of  the  69  squares  as  shown  in  Fig.  1.  are  included  in  the  region 
of  integration  for  the  corresponding  adjacent  squares  according  to^ 

H.  (2)  = J J H,(r,y)dxdy 
1 sJJa  1 

• • 

• • 

• ■*  • ,5 
H. (5)  s J 1 H (x,y)dxdy  + \ J J H (x,y)dxdy 

1 S5  1 D 

^(6)  ■ J j ^{x.yjdxdy  + | J I H^x.yjdxdy  (9) 

Sg  CUE 

and  so  forth  around  the  disk.  For  completeness,  the  underlying  basis  functions 
H^(x,y),  i«l,2,...,8  are  defined  in  Appendix  I . 


It  should  be  noted  that  the  basis  functions  H, (x,y),  i**l,2,...,8  vanish 

identically  outside  the  disk  of  radius  r»  *^69 pixels. 
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For  each  application  of  the  Hueckel  operator  an  edge  is  indicated  in  a 
4x4  block  by  marking  each  pixel  whose  center  is  within  ^2  pixels  of  the 
declared  edge.  The  computations  required  to  carry  this  out  have  not  been  in- 
cluded as  part  of  the  computation  count.  Since  the  operator  must  be  applied 
once  for  each  group  of  16  pixels  the  total  number  of  operations  for  each  appli- 
cation has  been  divided  by  l6  to  produce  the  computational  costs  per  pixel  as 
indicated  in  Table  4.  The  range  of  values  indicated  for  the  31  bit  arithmetic 
is  the  result  of  conditional  branches  within  the  program. 


Operation 

Required/Pixel 

16  bit  addition 

multiplication 

23.44 

6.94 

31  bit  addition 

multiplication 
square  root 

1.63-3.25  (average  2.44) 
2.5  -4.63  (average  3.57) 
.125-1. 56 (average  0.84) 

TABLE  4 

Computational  Costs  Per  Pixel  for  Hueckel  Operator 


Using  the  average  values  from  Table  4 together  with  the  normalized  computational 
costs  in  Table  1 we  arrive  at  an  effective  cost  of  109  operations  per  pixel. 

Two  parameters  CONF  and  DIFF  were  included  in  the  descriptions  provided 
by  Hueckel  [3]*  [4].  There  is  some  ambiguity  as  to  the  physical  significance  of 
these  parameters  but  generally  CONF  sets  the  desired  confidence  level  while  DIFF 
sets  the  smallest  level  difference  to  be  accepted  as  an  edge.  A program  listing 
describing  the  implementation  of  the  Hueckel  operator  employed  in  this  work  is 
provided  In  Appendix  III  from  which  the  exact  useage  of  these  parameters  can  be 
established. 
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V.  Results 


The  Wiener  and  Hueckel  operators  hare  been  exercised  on  two  types  of  images 


digitized  versions  of  actual  scenes  and  synthetic  images  with  additive  white 
Gaussian  noise  (i.e.,  independent  from  pixel-to-pixel ) . In  all  cases  we  have  in- 
cluded results  for  the  U-point  Laplacian  operator  strictly  for  comparison  purposes 


The  particular  images  considered  here  were  chosen  because  they  exhibit  performance 


that  is  typical  of  what  has  been  observed 


The  times  required  for  execution  of  the  two  programs  are  not  sufficiently 


different  to  be  noteworthy.  The  Hueckel  operator  was  divided  into  two  parts: 
projection  onto  the  basis  vectors  (execution  time  25  seconds)  and  computation 
(65  seconds)  for  a total  time  of  90  seconds.  The  Wiener  operator  likewise  has 
been  divided  into  two  parts:  linear  filtering  (56  seconds)  and  detection  (28 

seconds)  for  a total  time  of  dU  sc  "’•ds. 


The  performance  of  the  various  edge  detectors  is  Illustrated  in  Fig.’s  2 


through  7 vhen  the  input  Image  consists  of  the  previously  described  stochastic 


ed'e  process  imbedded  in  a white  Gaussian  noise  field.  In  all  cases  the  para- 
meters of  the  Wiener  operator  (T,F)  and  that  of  the  Hueckel  operator  (CONF,  DIFF) 


have  been  empirically  adjusted  to  optimize  performance.  The  Laplacian  has  been 


consistently  thresholded  on  the  positive  peak.  This  causes  some  displacement  of 


edges  toward  the  brighter  regions 


It  is  clear  from  Fig.'s  2 and  3 that  in  the  absence  of  noise  (i.e.,  S*00) 


both  the  Wiener  and  Laplacian  detectors  perform  well  when  there  is  a pronounced 
change  in  gray  level  across  an  edge  (i.e.,  large  negative  correlation)  such  as 


would  be  the  case  with  p*-0.9.  The  Hueckel  operator,  on  the  other  hand,  exhibits 


a curious  inability  to  resolve  closely  spaced  edges  and  a tendency  to  provide 


spurious  edge  declarations  particularly  in  regions  of  high  edge  density.  This  prob- 
lem becomes  more  pronounced  as  X,  the  average  edge  density /pixel,  increases. 


The  situation  is  further  aggravated  in  the  presence  of  less  pronounced  edges 


corresponding  to  an  Increased  correlation  p across  an  edge.  This  is  clearly 


illustrated  in  Fig.  k for  the  case  X«  0.05  and  p-0.5  still  with  C*».  It  should 
be  noted,  however,  that  the  Hueckel  operator  is  being  used  strictly  in  the 
edge  mode.  Had  it  been  used  in  the  line  or  edge-line  mode  it  is  possible 
that  this  problem. might  not  be  as  pronounced. 

The  behavior  as  a function  of  SNR  C with  p«-0.9  is  illustrated  in  Fig.»s  5 
and  6.  Here  X»0.0125  so  that  with  this  low  value  of  edge  density  the  edge 
resolution  questions  raised  above  are  not  an  issue.  Observe  that  for  £<10dB 
the  ordinary  Laplacian  is  virtually  useless  in  discerning  edge  structure. 

The  situation  becomes  even  worse  with  increasing  correlation  p as  illustrated 
in  Fig.  7 for  p«0.5.  The  Wiener  and  Hueckel  operators,  on  the  other  hand,  are 
both  able  to  detect  considerable  edge  structure  while  preserving  a considerable 
amount  of  noise  immunity.  Indeed,  the  Hueckel  operator  possesses  excellent 
immunity  to  this  type  of  salt-and- pepper  noise. 

Experiments  were  also  conducted  on  a synthetic  image  consisting  of  a light 
circle  (gray  level  6k)  intersecting  a polygonal  figure  (gray  level  31)  against 
a white  Gaussian  noise  field.  The  SNR  is  here  defined  in  terms  of  the  ratio  of 
gray  level  within  the  circle  to  the  standard  deviation  of  the  background  noise. 
Typical  results  are  illustrated  in  Fig.  8 for  SNR*6dB  and  selected  parameter 
choices  for  both  the  Wiener  and  Hueckel  operators.  Again  the  Laplacian  is 
virtually  useless  in  discerning  edge  structure.  Both  the  Wiener  and  Hueckel 
operators  exhibit  excellent  noise  immunity  with  the  Hueckel  operator  offering 
some  advantage  in  this  regard.  Furthermore  ,the  noise  immunity  of  the  Hueckel 
detector  is  observed  to  be  a definite  function  of  the  parameter  DIFF  by  comparing 
Fig. '8  8c. and  8d.  The  noise  immunity  properties  of  the  Wiener  detector  depend 
in  a more  complicated  way  upon  the  parameters  defining  this  operator.  Parameter 
choices  illustrated  in  Fig.  8 were  determined  empirically  to  be  near  optimum. 

Observe  that  the  Wiener  operator  has  a tendency  to  produce  connected  edge  segment 

1 

while  the  Hueckel  operator  tends  to  produce  disconnected  or  fragmented  edge  segments. 
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A Number  of  experiments  have  also  been  conducted  vith  various  classes 
of  real-world  images.  Results  for  a typical  head-and-shoulders  image  are 
illustrated  in  Fig.  9 while  results  for  a typical  chest  X-ray  are  shown  in 
Fig.  10.  Observe  the  flexibility  afforded  by  the  Wiener  detector  in  that 
by  logical  assignment  of  parameter  values  one  may  extract  different  portions 
of  the  "relevant”  edge  structure  in  typical  real-world  images.  This  is  demon- 
strated clearly  in  Fig.  10  and  to  a lesser  extent  in  Fig.  9*  In  particular, 
from  Fig.'s  lOe  and  lOf  note  that  the  Wiener  edge  detector  can  be  "timed"  to 
extract  the  closely  spaced  and  sharp  edge  structure (as  in  Fig.  lOe  by  the  para- 
meter choice  A*0.05  and  p«-0.9)or  the  less  frequent  and  more  gradual  edge  . 
structure (as  in  Fig.  lOf  by  the  parameter  choice  X*0.0125  and  p*0.5).  The 
former  would  be  useful  If  one  were  Interested  in  rib  structure  while  the  latter 
would  be  useful  if  instead  one  were  Interested  in  merely  extracting  the  outline 
of  the  chest  cavity.  Corresponding  results  are  illustrated  in  Fig.'s  10c  and  lOd 
for  the  Hueckel  operator  where  now  the  parameters  CONF  and  DIFF  have  been 
empirically  adjusted  to  emphasize  fine  and  coarse  structural  detail  respect- 
ively. Observe  the  inability  of  the  Hueckel  operator  to  discern  closely 
spaced  edges,  the  extreme  sensitivity  to  spurious  image  detail  and  the  tendency 
to  produce  disconnected  edge  segments.  The  Wiener  operator  is  far  superior  in 
this  regard.  Similar  consents  apply  to  Fig.  9> 

Finally,  in  Fig.'s  11  and  12  the  performance  of  these  two  edge  detectors 
is  illustrated  for  two  images  exhibiting  extensive  and  predominant  edge  structure. 
Observe  again  the  superior  behavior  of  the  Wiener  operator  vis- Si- vis  the  Hueckel 
operator  in  its  ability  to  be  "tuned”  to  extract  either  fine  or  coarse  structur- 


al detail. 
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In  summary  the  important  features  of  the  Wiener  operator  when  applied  to 
real-world  images  are:  relative  insensitivity  to  background  noise  and/or 


spurious  image  detail;  the  ability  to  discern  between  soft  edges  or  shading 
(gradual  transitions)  and  hard  edges  (step  transitions);  the  ability  to  dis- 


criminate between  closely  spaced  and  more  widely  spaced  edge  structure;  and  a 


tendency  to  produce  connected  edge  segments  which  segment  the  image.  These 
properties  have  been  shown  to  be  of  value  in  a contour  tracing  algorithm  [6] 


for  which  the  Wiener  detector  has  been  used  as  a front  end 


ator,  on  the  other  hand,  is  quite  sensitive  to  spurious  image  detail  and  suffers 
from  an  inability  to  discern  soft  edges  and/or  closely  spaced  edge  structure. 

A larger  disk  size  and/or  less  coarse  application  (i.e. , on  smaller  than  U x 
blocks)  might  help  but  only  at  the  expense  of  cm  increase  in  computational 


complexity . The  Hueckel  operator  also  tends  to  produce  disconnected  edge 


segments  thus  requiring  more  higher  level  processing  than  the  Wiener  detector's 


output  does  in  order  to  produce  a contour 


VI.  Conclusions 


utatlonal  complexity  and  suitability  for  specific  applications.  Problems  in 


making  these  comparisons  were  discussed  and  the  implementations  of  both  operators 

",  ■ ' : ' 

used  for  this  comparison  have  been  described  in  detail. 


Both  real  world  and  synthetic  Images  were  considered  for  the  purpose  of 


evaluating  performance.  While  the  Hueckel  operator  exhibited  better  subjective 
imunity  to  salt-and-pepper  noise  than  the  Wiener,  the  Wiener  operator  excelled 
in  its  ability  to  discern  between  shading  and  edges.  In  addition,  the  Wiener 
operator's  tendency  to  produce  connected  boundaries  make  it  more  suitable  as 
a front-end  processor  for  a contour  tracing  algorithm. 


Appendix  I 

Basis  Functions  for  Hueckel  Operator 


The  basis  functions  H^x.y) , 1*1,2, . . . ,8  from  which  the  discretized  versions 
H^j),  i*l»2,...,8,  J*l,2,... ,69  are  computed  according  to  (8)— (9)  are  defined 

as  follows.  Let  (x,  y)  represent  the  position  relative  to  the  center  of  the 
disk  in  pixels.  The  disk  is  assumed  of  radius  r * /*Vir  pixels.  Define 
normalized  coordinates  x * A/m  x and  similarly  y * A/tt  y. 


52  * x2  + y2 
and  , , 

Set  H^x.yJ-O  forfi2  >1,  i*l,2,...,8  while  for  62  < 1: 


(1-1) 

(1-2) 


Hi(x»y) 

H2(x,y) 
Hg(x,y ) 
^(x.y) 
H5(x,y) 
Hg(x,y) 
H7(x,y) 


Q(5«2-2) 


AS  Q(x2-y2) 
AS  Q(2xy) 
AS  Q(262-1)x 
AS  Q(262-l)y 


H8(x,y)  • Q( -2+176 2-2l6") 


(1-3) 
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AVAILABLE  COPY 

Appendix  II  — • 

Program  to  Implement  Wiener  Operator 


NPIX  IS  THE  SIZE  OF  THE  IMAGE 
X IS  THE  INPUT  IMAGE 
V IS  THE  FILTERED  IMAGE 

E IS  AN  ARRAY  OF  FLAGS:  0 FOR  NO  EDGE;  -1  FOR  AN  EDGE 

YL  AND  YR  ARE  ARRAYS  USED  TO  HOLD  CURRENT  AND  PREVIOUS 
OUTPUTS  OF  THE  LEFT  AND  RIGHT  MOVING  FILTER  SECTIONS. 

DIMENSION  X<NPIX,  NPIX>,  Y<NPIX,  NPIX>,  E<NPIX, NPIX> 

DIMENSION  VL<2,  NPIX>,  YR<2, NPIX) 

CUR  AND  PREV  POINT  TO  THE  CURRENT  AND  PREVIOUS  PARTS  OF  YL  AND  VR 
THEY  ARE  INITIALIZED  BELOW. 

INTEGER  E,  CUR,  PREV 

CUR-1 

PREV-2 

Bll, A10,  All,  AND  A ARE  PARAMETERS 
OF  THE  FILTER  THAT  MUST  BE  ENTERED 

F AND  TRESH  ARE  PARAMETERS  OF  THE  EDGE  DETECTION  ALGORITHM 
READ  Bll>  A10,  All,  A,  F,  THRESH 
BIG— . 5*<Bil*l.  ) 

INITIALIZE  PREVIOUS  LINE  ARRAYS  TO  ZERO 
DO  100  J— 1, NPIX 
VL<PREV,  J)-©. 

YL<CUR,  J >— 0. 

YP/.PREV,  J>-0. 

YR<CUR,  J>»0. 

CONTINUE 

SET  FIRST  LINE  OF  FILTERED  IMAGE  TO  © 

DO  200  J— 1,  NPIX 
Y<1,  J>«0. 

CONTINUE 

MOVE  DOWN  THROUGH  IMAGE  LINE  BY  LINE 
DO  60©  1-2, NPIX-1 

MAKE  OLD  CURRENT  LINE  NEW  PREVIOUS  LINE 
I HOLD-CUR 
CUfi-FREV 
PREV- I HOLD 

SET  FIRST  AND  LAST  POINTS  O*’  FILTERED  IMAGE  TO  0 
Y<  I,  l>-0. 

V<  I,  NPlX>-0. 

MOVE  THROUGH  THE  LINE  POINT  BY  POINT 
DO  40©  J-2,  NPIX-1 


IT 


* 


II 


BEST  AVA1UBLE  COPY 


C RIGHT  MOVING  FILTER  SECTION 

YR'CUR-  J>=-A10*<YR<CUR/  J-1>+YR<PREV,  J> >-All*YR<PREV,  J-l> 

C vx<I, J>+B10*<<X<<I/ J-1>+X<I-1,  J>>+B11*X<I-1,  J-l> 

C LEFT  MOVING  FILTER  SECTION 

JL-NPIX-J+1 

VL<CUR,  JL>— A10*^L<CUR/  JL+1>+YL<PREV,  JL>>A11*YL<PREV,  JL+1> 

C +XCI,  JL>+B10*<X<I,  JL+l>Vx<I-l,  JL>>+B11*X<1-1/  JL+1> 

C END  POINT  BV  POINT  LOOP 

400  CONTINUE 

C SUM  LEFT  AND  RIGHT  MOVING  FILTER  OUTPUTS 

DO  500  J-2,  NPIX-1 
V<I,  J>«YR<CUR,  J>+YL<CURj J+l> 

000  CONTINUE 

C END  LINE  BV  LINE  LOOP  MOVING  DOWN 

600  CONTINUE 

C ZERO  LAST  LINE  OF  FILTERED  OUTPUT 

DC*  706  J«l,  NPIX 
YCNPIX,  J>»0. 

700  CONTINUE 

C GET  READY  TO  MOVE  UP  THROUGH  THE  IMAGE 

C ZERO  PREVIOUS  LINE  ARRAY 

DO  11O0  J»l>  NPIX 
VL<CUR,  J>«0. 

YL'XPREV,  J>«0. 

VP. < CUR/  J>®0. 

VP.  < PREY/  J>»0. 

1100  CONTINUE 

C MOVE  UP  THROUGH  THE  IMAGE  LINE  BY  LINE 

DO  1600  II-2/NPIX-2 
I-NPIX+l-II 

C MAKE  OLD  CURRENT  LINE  NEW  PREVIOUS  LINE 

I HOLD-CUR 
CJR-PREV 
PREVMHOLD 

C MOVE  THROUGH  THE  LINE  POINT  BY  POINT 

DO  1400  J«2/ NPIX-1 

C RIGHT  MOVING  FILTER  SECTION 

YR<CUR>  J>«-A10*<YR<CUR,  J-1>+VR<PREV/  J>>-A11*YR<PREV/  J-l> 

C *X<I/  J>+B10+<<X<<I/  J-1>*X<I*1.  J>>*B11*X<H-1,  J-l> 

C LEFT  MOVING  FILTER  SECTION 


■J 


■aem* 


JL-NPIX-J+1 

VL<CUR,  JL>— A10*(^L<CUR,  JL+1>«-VL<PREV,  JL>>A11*YL<PREV,  JL+1> 

C +X<I,  JL>*B10*<X<I,  JL*1>+X<  1+1/  JL>>+B11*X<I+1,  JL+1> 

END  POINT  BV  POINT  LOOP 
CONTINUE 

ADD  OUTPUT  OF  LEFT  DOWN  PLUS  RIGHT  DOWN  TO  LEFT  UP  PLUS  RIGHT  UP 
AND  MULTIPLY  SUM  BV  A 
CO  15B0  J“2/ NPIX-1 

V<  1-1/  .t)*A*(Y<I"L  J>+VR<CUR,  J>+YL<XUR,  J+l>> 

CONTINUE 

END  LINE  BV  LINE  LOOP  MOVING  UP 
CONTINUE 

END  OF  FILTERING  - START  OF  DETECTION 

2ER0  FIRST  AND  LAST  LINES  OF  EDGE  DETECTOR  OUTPUT 

DO  2100  J*l,  NPIX 

E<I,  J>«0 

E<NPIX,  J>«0 

CONTINUE 

MOVE  THROUGH  THE  IMAGE  LINE  BV  LINE 
DO  2400  1*2, NPIX-1 

ZERO  EDGE  DETECTOR  OUTPUT  FOR  FIRST  AND  LAST  POINTS 

E<I,1>“0 

E<X#NPlX>*e 

MOVE  THROUGH  THE  LINE  POINT  BV  POINT 
DO  2500  J“2, NPIX-1 
D I FI* ABS < V < 1-1,  J-i>-v<<I+l, J+l>  > 

D1F2“ABS<V<1-1,  J>-V<I«-1,  414 

DIF3*RBS<V<  I-i, J+l>-V<  1+1,  J-l>  > 

DIF4»ABS<V<<I,  J-l>-V< It J+1>>*1.  414 
DIF*AMAX<OIFl,  DIF2,  DIF3/DIF4> 

SUM*V< I,  J-1>*V<I,  J)+V<L  J+l> 

C *.V<I-1,  J-1>+V<I-1,  J>+V<I-1,  J«-l> 

C ♦v< i+i,  jr-i>+v< i+l,  J >+v<  i+i,  J+l> 

SUM«SUM*SUM*F 

FLAG*© 

IF  <DIF**2-SUM.  GT.  THRESH > FLAG— 1 
E<ii  J>«FLAG 

END  POIHT  BV  POINT  LOOP 
CONTINUE 

END  LINE  BV  LINE  LOOP 

CONTINUE 

STOP 
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Appendix  III 


Program  to  Implement  Hueckel  Operator 
A)  16-Bit  Arithmetic  Used  in  Calculating  A1-A8 


♦HK  1>*<  X<  1>*X<12>*X<49>*X<65>*X<69>*,X<57>+X<21>*X<  5>  > 
+H1<  2>*<  X<  2>*X<20>*X<48>+X<$6>+X<€8>«-X<4e>«-X<22>«-X<  4>  > 
♦Hl<  3>*<  X<  3>«-X<31>*X<67>-*-X<39>  > 

♦MIC  «>*<  X<  6>^X<12>+X<64>+X<58>  > 

>H1<  7>*<  X<  7>*X<ll>+X<14>*X<50>+X<63>«-X<59>«-X<56>*X<r:20>  > 
♦Hl<  8>*<  X<  8>+X<10>+X<29>^X<47>+X<62>^X<80)*X<41>^X<22)  > 
«-Hi<  ?>*<  X<  9>-*-X<32>+X<61>+X<28>  > 

X<19>^X<15>^X<51>^X<55>  > 

X<18>*X<16>+X<28>+X<46>+X<52>*X<54>«-X<42:>«-X<24>  > 
♦Hl<  17>*''.  X<i.7>*X<33>*X<53>*X<3?>  > 

♦Hl<25>*r  X<23>*X<27>*X<45>+X<43>  > 

♦Hl<26>*<  X<26>-*-X<34>-»-X<44>*X<3€>  > 

•*,Hi.<'.35>»*<<  X<35>  > 

♦H2<  5>*<  X<  5>*X<69>-X<  1>-X<65>  > 

♦H2<21>*<  X<21>-*X<57>-X<13>-X<49>  > 

«-H2<  4>*<  X<  4>«-X<68>-X<66>-X<  2>  > 

♦H2<22>*<  X<22>-*X<40>-X<20>-X<48>  > 

♦H2<3S*>f<  X<39>-X<31>  > 

♦H2<  €>*<  X<  8>*X<38>-X<12>-X<64>  > Rj-\T  Ay  All  AR 

■*H2<  7>*<  X<  7>-*X<59>-X<ll>-X<63>  > ULJ I _ M ¥ r\l  LnUI 

«-H2<19>*<  X<19>^X<55>-X<15>-X<51>  > 

♦H2<20>*<  X<20>^X<56>-X<14>-X<50>  > 

■*-H2<  8>*<  X<  8>-*-X<60>-X<10>-X<62>  > 

♦H2<23>*<  X<22>*X<41>-X<29>-X<47>  > 

•*H2<28>*<  X<38>-X<32>  > 

♦H2<18>*<  X<18>+X<54>-X<16>-X<52>  > 

♦H2<24>*<  X<24>+X<42>-X<26>~X<46>  > 

♦H2<37>*<  X<37>-X<33>  > 

♦H2<25>*<  X<25>*X<43>-X<27>-X<45>  > 

«-H2<36>*<  X<3€>-X<34>  > 

+H3<  1>*<  X<  i>«-X<  5>“X<65>“X^89>  > 

+H3<21>*<  X<21>+X<13>-X<49>-X<57>  > 

♦H3<  2>*C  X<  2>+X<  4>-X<66>-X<68>  > 

♦H3<22>*<  X<22>-t-X<30>-XC48>-X<40)  > 

♦H3<  3>*<  X<  3>-X<67>  > 

♦H3<  6>*<  X<  € > ♦X < 12 > -X < 64 > —X < 58  > > 

♦H3<  7>*<  X<  7>**-X<ll>-X<63>-X<59>  > 

♦H3<20>*<  X<20>+X<14>-X<56>-X<50>  > 

♦H3<  8>*<  X<  8>*t-X<10>-X<62>-X<€0>  > 

♦H3<23>*<  X<23>«-X<29>-X<4?>-X<41>  > 


+H3<  9>*<  X<  9>-X<61>  > 

•*-H3<19>*<  X<19>+X<±5>--X<51>*-X<35> 
■*-H2<:i8>*<  X<18>+X<16>-X<52>-X<54> 
♦H3<24>*<  X<24>+X<28>-X<46>-X<42> 
«-H3C17>*'  XC17>-X<52>  > 

♦H3<25>*<  X<25>+X<27>-X<45>-X<43> 
«-H2<26:>*‘’.  X<26>-X<44>  > 

■*-H4'.21>*<  X<21>«-X<13>+X<49>-*-X<57>- 
«-H4<22>*<  X<22>+X<38>-t-X<48>^X<40>- 
+H4<39>*<  X<39>+X<31>-X<67>-X<  3> 
+H4<20>*<  X<20?-»-X<14>-*-X<50>+X<5€>- 
+H4C23>*<  X<23>+X<29>+X<47>+X<41>- 
-»-H4<38>*<  X<38>+X<32>-X<61>-X<  9> 
♦H4<24>*<  X<24>-*-X<28>-t-X<46>+X<42>- 
♦H4<37>*<  X<3?>*X<33>“X<53>-X<17> 
+H4<36>*<  X<36>+X<34>-X<44>-X<26> 
♦H5<  5>*<  X<  5>+X<49>*X<65>*X<21>- 
♦HSC  4>*<  X<  4>+X<48>-t-X<66>-t-X<22>* 
■*‘H’3<  «>*<  X<  6>+X<64>-X<12>-X<58> 
♦H3<  7>*<  X<  7>*X<63>+X<50>*X<20>- 
♦H5<  8>*<  X<  8>+X<47>+X<€2>+X<23>- 
♦H5<19>*<  X<19>*X<51>-X<15>-X<55> 
*H5<18>*<-  X<18}+X<46}*X<52>+X<24>* 
♦H5'.25>*<  X<25>+X<4,5>-X<27>-XC43> 
+H6C  5>*<  X<  5>*X<69>-X<  1>-X<C65> 
■*-HS<21>*<  X<21>«-X<37>-X<13>-X<49> 
*+k*<  4>*<  X<  4>*X<68>-X<  2>-X<66> 
■*+46' 22>*<  X<22>^X<4O>-X<30>-X<48> 
♦H6<39>*<  X<39>-X<31>  > 

■*■«€<  6>*<  X<  6>*X<58>-X<12>-X<64> 
♦H6<  7>*<  X<  7>*X<59>""X<11>“X<63> 
■*-H6<19>*<  X<19>+X<55>-X<51>-X<15> 
■*;H€<20>*,<  X<20>+X<56>-X<14>-X<50> 
+H€<  8>*<  X<  8>+X<60>-X<10>-X«S2> 
+Hf'.a3>*<  X<23>+X<41>-X<29>-X<47> 
■»+l6<38>*<  X<38>-X<32>  > 

♦H6<18>*<  X<18>-*-X<34>-X<16>-X<52> 
+H€<24>*<  X<24>*X<42>-X<28>-X<48> 
♦H€<37>*<  X<37)-X<33>  > 

♦H6<25>*<  X<23>-0«43>-X<27>-X<45> 


X<65> 

X<€6> 

> 

X<11> 

X<10> 

> 

XC16> 

> 

> 

X<  1> 
X<30> 
> 

•X<11  > 
•X<10> 
> 

■X<16> 

> 

> 

> 

> 


-XC69>— X<  5>- 
-X<68>-X<  4>- 

-XC€3>-X<59>- 

-X<€2>-X<60>- 


X<  1>  > 
X<  2>  > 

•X<  7>  > 
•X<  8>  > 


-X<52>-X<54>-X<18>  > 


■X<13>-X<€9> 

■X<68>-X<40>- 

•X<14>“X<59>" 

-X<29>-X<60>- 


XC57>  > 
X<  2>  > 

■X<56>  > 
•X<41>  > 


-X<28>-X<54>-X<42>  > 
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+H€<36>*<  X<36>-X<34>  > 
ft?  « +H7<  1>*<  X<.  1>+X<  5 > -X  < 65  > —X  < 69  > > 

+H7<21>*<  X<21>+X<13>-X<49>-X<57>  > 

♦H7<  2>*<  X<  2>+X<  4>-X<66>-X<68>  > 

+H7<22>*<  X<22>+X<30>~X<48>”X<40>  > 

+H7<  3>*<  X<  3>-X<67>  > 

+H7'  €>*'.  X<  6>+x<12>-X<64>-X<58>  > 

+H7<  7>*<  X<  7>+X<ll>-X<63>-X<59>  > 

r*H7<20>*<  X<20>+X<14>-X<56>-X<50>  > 

♦H7<  8>*<  X<  8>*X<10>-X<62>-X<60>  > 

+H7<23>*<  X<23>+X<29>-X<47>-X<41>  > 

+H?<  9>*<  X<  9>-X<61>  > 

♦H7<19>*<  X<19>+X<15>-X<51>-X<55>  > 

«-H?<i8>*<:  X<18>*-X<±6>-X<52>-X<34>  > 

♦H7':24>*<  X<24>+X<28>-X<46>-X<42>  > 

+H7<17>*<  X<17>-X<53>  > 

♦H7<25>*<  X<25>+X<27>-X<45>-X<43>  > 

♦H7< 26 >*<  X<26>-X<44>  > 

fi8  - -t-H8<  1>*<  X<  1>+X<13>«-X<49>+X<65>+X<69>+X<57>+X<21>+X<  5>  > 
♦H8<  2>*<  X<  2)+X<30>+X<48>+X<66>+X<68>+X<40>*X<22>*X<  4>  > 
+H8<  3>*<  X<  3>+X<31>+X<67>+X<39>  > 

•‘■H8''.  6>*<  X<  6>+X<12>+X<64>+X<58>  > 

«*H8<  7>*<  X<  7>+X<ll>-*-X<14>+X<50>+X<63>+X<59>-*-X<56>*X<20>  > 
+H8<  8>*<  X<  8>+X<10>+X<29>+X<47>*X<62>*X<60>+X<41>+X<23>  > 
+H8<  ?>*<  X<  9>+X<32>+X<61>+X<38>  > 

+H8<19>*<  X<19>+X<15>+X<51>+X<55>  > 

♦H8<18>*<  X<18>+X<16>+X<28>+X<46>+X<52>+X<54>+X<42>+X<24>  > 
♦H8<17>*<'.  X<17>+X<33>+X<53>*X<37>  > 

+H8<20>*<  X<23>+X<27>+X<45>+X<43>  > 

+H8<26>*<  X<26>+X<34>+X<44>+X<36>  > 

+H8<35>*<'.  X<35>  >. 


B£S1  available^® 
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B)  31-Bit  Arithmetic  Used  in  Fitting  Ideal  Edge  Element 


REfiL  ftl.  ft2.  ft3.  ft4.  ft5.  ft6.  ft?.  ft8.  ft2€..  ft27.  ft27.  IPU2 
PERL  El..  E2-  E3.  E4.  E5.  E6.  Gl-  G2,  G3.  CX.  OCX.  CV,  CCV 
REftL  MU.  NU.  00.  01.  U2.  U3.  Kft.  W.  PP.  BP..  DP.  DIFF.  DIFFS 
PERL  ED  JS.  TP.  CONF.  CONFM1  - RM.  BM.  DM.  TM.  THRESH.  MflXBRT 
THRESH  = . 0001 
MftXBRT  *»  25. 

DIFFS=DIFF**2 
C0NFMi*l-CONF 
fl26*ft2**>2+ft6**2 
R37«ft3**2+ft?**2 
ft2?»fl37+ft26 

lP02=ftl**2+ft27+ft4**2+ft5**2+R8**2 
IF  IPU2.  LT.  27*PI/64  GO  TO  DISft 
El«<l/3>**.  5*ftl 
E2«<:2/3>**.  5*ft4 
E3“<2/3>**.  5*fi5 
E4“2*E1*E2+ < ft26-ft2?  > /2 
E5*=2* El*E3+ft2*fi3+fl6*fi7 
G1“SIGN< E2*E4+E3*E5  > *SQRT  < E2**2+E3**2  > 

CCX*G1*E2+E4 
CCV^G1*E3+E5 
G2-SQRT  ‘:CCX**2-*-CCV**2> 

CCX*CCX/G2 

CCV-CCV/G2 

E6*E2*CCX+E2*CCV 

SIG2-El**2+ft27/2+fl8**2+E6**2+E4*CCX+E5*CCY 

IF  SIG2/IPU2-K1-C0NF>*SIG2**2/<SIG2**2+DIFF**2>.  LE  1 0 GO  TO  DISft 
G3"<2*<1+CCX> >**.  5 

IF  G3-THRESH<0  THEN  CV*=SIGN<CCV>  ELSE  CY*CCY/G3 
CX-G2/2 

U0»ft2*CX+ft2*CY 
01“ < El+ES >*3**-l/2 
U2^  <.  ft€*CX+ft?*CV  > *5**-.  3 
U2*  •:  Ul-2*ft8  > /5 
NU»U1**2-U2*U0 

IF  ftBS<  NO) -THRESH,  GT.  0 GO  TO  NTRIC 
IF  flBS<MO>.  LE.  THRESH  GO  TO  DISft 
IF  ftBS<U0>.  GE.  THRESH  GO  TO  TftU0NZ 
IF  ftBS<01>.  GE  THRESH  GO  TO  TftUINZ 
IF  ftBS<02>.  LT.  THRESH  GO  TO  DISft 


BEST  AVAILABLE  COPT 


t 


bl  AVAILABLE  CORY 


KA-U2VU2 


GO  TO  KTEST 
TRO0N2  KR-Ul/UO 

GO  TO  KTEST 
T RUINS  KA-U2/U1 

GO  TO  KTEST 
N7RIC  H*J-U2**2-U1*U2 

KR- < 01*02-00*02 > /2*N0 
H-Kfl**2-MU/NU 

IF  RBS  < W > - THRESH.  LT.  O 00  TO  KTEST 
IF  W.  LT.  O 00  TO  COMPW 
H»SQP.T<M> 

RM-KA-W 

RF-KR*W 

BM«1-RM**2 

BP«i-P.P**2 

IF  ABS<RM>.  QE.  . 87  . OR.  RBS<RP>.  QE.  . 87  00  TO  RLTRP 

0M-2<  2#  >**-.  5*  <O0*RP-O1 

TM«WVW 

0P-2< 2#  > **-.  5* < Oi-RM*O0 > /BP* *2 
TP»OP/W 

IF  RBS-:TM>.  QE  MRXBRT  OR  ABS^TPX  QE  MRXBRT  00  TO  RLTRP 
00  TO  TTEST 
RLTRP  TM-0 

RP-RM-Ui/00 

IF  flBS'!RP>.  QE.  0.  S»2  00  TO  OISR 
TP*U0/<  <1-RP**2>**2*SQ3P4> 

IF  ABS<TP>.  QE  MRXBRT  00  TO  OISR 
GO  TO  TTEST 

KTEST  IF  KR.  QE.  1 00  TO  01  SR 

RM-RP-Kfi 

TM*TP«U0*2/SQRT < 3*P I > / < 1-KR**2> **2 
TTEST  EOJS-TM+TP 

IF  ABS<TM*TP>.  QT.  MRXBRT  00  TO  RLTRP 

F:»  «C  RBS  < TM  > *RM*AB£  < TP  > *PP  > X f.  RBS  < TM  >+ABS<TP>  > 

IF  EOJS<>©  00  TO  CRIT 
EOJS*-- EOJS 

cx*— cx 

CY*-CY 
R— R 


RLTRP 


KTEST 


TTEST 


n 


CP. IT  GO  TO  EXIT 
COMF'H  R-KR 

EDJS-0 


PET'JPN 

FLRG-0 

RETURN 
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c)  Wiener  Operator  d)Hueckel  Operator 


Fig.  3 

Performance  of  Edge  Detectors  on  Random  Field  for  X=.05,  P=-.9,  C5”0’. 

(a)  Original  image;  (b)  Laplacian;  (c)  Wiener  operator  for  F=.3  and 
T=125;  (d)  Hueckel  operator  for  CONF».75  and  DIFF=1. 
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a)  Original  Image 


b)Laplacian 


c)  Wiener  Operator 


d)  Hueckel  Operator 


Fig.  h 

Performance  of  Edge  Detectors  on  Random  Field  for  X*.05,  P=.5,  ?=0°: 
(a)  Original  image;  (b)  Laplacian;  (c)  Wiener  operator  for  F=.9  and 
T*3«13;  (d)  Hueckel  operator  for  CONF*.75  and  DIFF»  1. 
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c)  Wiener  Operator  d)Hueckel  Operator 

n 

Fig.  5 
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Performance  of  Edge  Detection  on  Random  Field  for  X=,0125,  p=-.9,  X=10dB: 
(a)  Original  image;  (b)  Laplacian;  (c)  Wiener  operator  for  F=.5  and 
T»  (d)  Hueckel  operator  for  CONF-.95  and  DIFF=  9- 
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c)  Wiener  Operator 


d)  Haeckel  Operator 


Fig.  6 

Performance  of  Edge  Detectors  on  Random  Field  for  A=.0125,  p=-.9,  C=3dB: 
(a)  Original  image;  (b)  Laplacian;  (c)  Wiener  operator  for  F=,25  and 
T*  235.  (d)  Hueckel  operator  for  CORF* . 95  and  DIFF*  9. 


a)  Original  Image 


b)  Laplacian 


d)  Hueckel  Operator 


c)  Wiener  Operator 


Performance  of  Edge  Detection  on  Random  Field  for  A*. 0125,  p=.5>  C=3dB 
(a)  Original  image;  (b)  Laplacian;  (c)  Wiener  operator  for  F=.7  and 
T»  9.U;  (d)  Hueckel  operator  for  C01IF».T5  and  DIFF=  5. 


•)  Wi«fl«r  \ X-0.025,  p* O.j 
C*3,F«.I  ,T  *4000 


f ) Wl*n«r  \ X ■ O.OI25,  i*. 5 
{■3,F«.2,T»250 


Fig.  8 

Performance  of  Edge  Detectors  on  Synthetic  Image:  (a)  Original  image; 

(b)  Laplacian;  (c)  Hueckel  operator  for  CONF*.75  and  DIFF=  1;  (d)  Hueckel 
operator  for  CONF*.75  and  DIFF*  5;  (e)  Wiener  operator  for  A=.025,  p*0, 
C*3dB,  F».l  and  T*625;  (f)  Wiener  operator  for  X*.0125,  p=.5»  Cs3dB, 

F-.2  and  T«  3.9l. 
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c)  Hueckel  - Fin*  Detail  d)  Hueckel  - Coar**  Detail 


*)  Wiener  - Fine  Detail  f ) Wiener  - Coaree  Detail 


Fig.  9 

Performance  of  Edge  Detectors  on  Typical  Head  and  Shoulder  Image: 

(a)  Original  image;  (b)  Laplacian;  (c)  Hueckel  operator  for  CONF=.75 
and  DIFF«  1;  (d)  Hueckel  operator  for  COHF*.85  and  DIFF»  6;  (e)  Wiener 
-operator  for  A».025»  p«0,  C«3dB,  F*.5  and  T*l.'56;  (f)  Wiener  operator 
for  A*. 0125,  p*.5»  C*3dB,  F«.2  and  T«.08. 
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a) Original  Image 


b)  Lapiacian 


c)  Hueckel  - Fin*  Detail  d)  Hueckel  - Coarse  Dstail 


•)  Wiener  - Fine  Detail  f)  Wiener  - Coarse  Detoll 

Fig.  10 


Performance  of  Edge  Detectors  on  Typical  Chest  X-ray:  (a)  Original  image; 

(b)  Lapiacian;  (c)  Hueckel  operator  for  CONF*.75  and  DIFF=  1;  (d)  Hueckel 
operator  for  CONF*.85  and  DIFF»  6;  (e)  Wiener  operator  for  A-.05,  P*--9, 
(*10dB,  P^.5  and  T»1.56;  (f)  Wiener  operator  for  A«.0125,  p*.5,  C*3dB, 
F*.2  and  T».08. 
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a)  Original  Image 


b)  Loplaclan 


c)  Hueckel  - Fin*  Detail 


d)Hu*ckel-  Coarse  Detail 


a)  Wiener  - Fine  Detail  f) Wiener  - Cooree  Detail 

Fig.  11 

Performance  of  Edge  Detectors  on  Outdoor  Scene  1:  (a)  Original; 

(b)  Laplacian;  (c)  Hueckel  operator  for  CONF*.75  and  DIFF=  1;  (d)  Hueckel 
operator  for  CONF*.85  and  DIFF*6;  (e)  Wiener  operator  for  A=.05,  p=-.9, 
C»10dB,  F«.3  and  T*  125;  (f)  Wiener  operator  for  X«*.0125,  p=.5,  C=3dB, 
F«.2  and  T-  3.91. 
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a) Original  Image  b)  Laplaeian 


c)  Hueckel  - Fine  D#tail  d)  Hu#ck#l  - Coors#  Detail 


•)  Wl#n#r  - Fin#  Detail  f ) Wi#n#r  - Coart#  Detail 


Fig.  12 


Performance  of  Edge  Detection  on  Outdoor  Scene  2:  (2)  Original; 

(b)  Laplaeian;  (c)  Hueckel  operator  for  C0NF=.75  and  DIFF=  1; 

(d)  Hueckel  operator  for  CONF».85  and  DIFF*  6;  (e)  Wiener  operator 
for  A«.05,  p«-.9»  C*10dB,  F».3  and  T»  125;  (f)  Wiener  operator  for 
X-.0125,  p».5,  ;»3dB,  F-.2  and  T-3.91. 
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Image  Processing 
Edge  Detection 
Two-Dimensional  Filtering 


(TRACT  f< 


^*This  paper  is  concerned  with  an  empirical  evaluation  of  the  relative  per- 
formance of  several  selected  edge  detectors  for  the  detection  of  edge  structure 
in  noisy  digitized  images.  In  particular,  we  consider  both  the  Hueckel  operatcn 
and  a recently  introduced  class  of  edge  detectors  implemented  as  two-dimen siona: 
infinite  inpulse  response  (IIR)  or  recursive  digital  filters.  The  latter  were 


infinite  inpulse  response  (IIR)  or  recursive  digital  filters 
originally  developed  on  the  basis  of  least  mean-square  Wiener  spatial  filtering 
concepts  for  an  assumed  stochastic  model  of  edge  structure  in  typical  imagery  — 
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data.  These  two  schemes  are  compared  with  respect  to  performance,  computa- 
tional complexity  and  suitability  for  specific  applications.  Problems  in 
making  these  comparisons  are  discussed.  The  implementations  of  both 
operators  used  for  this  comparison  are  described  in  detail. 
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